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We study the static potential between external quark-antiquark pairs in a strongly coupled gauge 
theory with a large number of colors and massive dynamical flavors, using a dual string description. 
When the constituent mass of the dynamical quarks is set below a certain critical value, we find a 
first order phase transition between a linear and a Coulomb-like regime. Above the critical mass the 
two phases are smoothly connected. We also study the dependence on the theory parameters of the 
quark-antiquark separation at which the static configuration decays into specific static-dynamical 
mesons. 
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1. Introduction.- The study of non-perturbative ef- 
fects of dynamical quarks is a notoriously difficult prob- 
lem to address with the present computational techniques 
in QCD. Phenomenological models, such as the screened 
version [1] of the Cornell potential for heavy quarks [2], 
are of fundamental importance but an understanding of 
the physics of these effects from first principles is still 
missing. In fact, the best tool for the study of QCD at 
strong coupling, i.e. its lattice formulation, is still par- 
tially limited by the lack of suitable computational power 
when dealing with dynamical light flavors (for our pur- 
poses the relevant reference is [3]). In the latest years, 
string theory has developed as a possible tool for study- 
ing strong coupling effects in quantum field theory. It 
is clearly not an ideal setting, since it allows to study 
mainly the large 't Hooft coupling, large N c (number of 
colors) limit of theories which typically include a number 
of other fields besides (possibly) the QCD ones. Never- 
theless, apart from its obvious theoretical interest, string 
theory has proven to be a valuable way for understand- 
ing processes which are universal enough, such as specific 
Quark-Gluon Plasma ones. 

In this letter we study the "connected part" (i.e. we 
do not include the mixing with the meson/anti-meson 
states) of the static potential between two heavy test 
quarks in a SQCD-like quantum field theory, which in- 
cludes light dynamical quarks of mass m q , by using its 
dual string theory description. Below a critical mass m c 
of the dynamical quarks, we observe a first order phase 
transition in the potential V(L) as we increase the dis- 
tance L between the two test quarks. The transition 
takes place in the region where the potential turns from 
Coulomb-like to linear. Instead, for masses larger than 
m c , the connection between the two regions is perfectly 
smooth. We also study the "string breaking length" L s {, 
at which the heavy quark configuration decays into a spe- 
cific pair of heavy-light mesons. It is shown that L s b is 
a decreasing function of the number Nf and the mass 



m q of the dynamical flavors, the latter being a fully non- 
perturbative effect. 

While the theory we study is not QCD, it is possi- 
ble that the main features of the static potential are 
sufficiently universal. Hence, we provide an analysis of 
non-perturbative effects of light dynamical quarks in a 
strongly coupled quantum field theory, in a very straight- 
forward and simple way. 

Sections 2 and 3 provide the technical details on the 
string set-up. Section 4 and 5, which can be read inde- 
pendently, contain the results. We define x = jf- and 
limit our analysis to the range < x < 1.5 where we 
reach sufficient numerical precision. 

2. Details on the dual string description.- The 

SQCD-like theory we consider is realized at the four 
dimensional intersection of N c 3> 1 "color" D5-branes 
wrapped on a S 2 and Nf ~ iV c "flavor" D5-branes. It is 
a particular version of M — 1 SQCD deformed by a quar- 
tic superpotential as considered in [4]. In these models, 
the flavor branes are homogeneously smeared along the 
(large) internal manifold [5] . The typical separation be- 
tween two branes is much larger than the string scale so 
the flavor symmetry is in fact U(l) Nf . We find solutions 
of IIB supergravity coupled to the flavor brane sources. 
These sources are represented as a sum of DBI+WZ ac- 
tions [6] . Notice that even if the DBI action for Nf branes 
is not valid when g s Nf is not small, here we can use it 
because the branes are separated far apart so what we 
have is a sum of actions for individual branes. The string 
frame metric reads 
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The dilaton <f> and the functions Y, H, G, a depend only on 
the radial coordinate p; we used the notation g s = e 1 ^" 1 , 
where (pin = (p(p = 0). In the massless flavor case, which 
was studied in [4] , the Nf smeared branes were extended 
along p from the origin up to infinity. In the following 
we will consider the case in which all the flavor branes 
are extended up to a finite distance p q > from the 
origin. This will correspond to turning on a mass for all 
the flavors in the gauge theory dual. 

The full supergravity background contains also a RR 
three-form given by 
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where a> 3 = dtp + cos 9 'dip + cos 9 'dip. In [4], a set of dif- 
ferential equations and constraints for H, G, b, a, Y, (p was 
found. They read, in the present notation, 

TT C(a + b) 2-x „ C(a-b) 2-x 

H = , G = 1 , 

4 8 ' 4 8' 
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where C = cosh(2p), S = sinh(2p). This system of equa- 
tions is valid for arbitrary Nf, N c . 

The expressions (2) make the addition of massive fla- 
vors relatively easy. As it was noticed in a similar context 
[7], this can be effectively achieved by just replacing the 
constant x with a function of p, x(p). In this way the 
expression for the three-form flux is unchanged, modulo 
the above redefinition of x. Indeed the Bianchi identity 
gets changed, but since in the supersymmetric variations 
for the background fermions only the field F (and not its 
derivative dF) enters, the whole BPS equations will not 
be modified in form. On the field theory side, we can 
neglect the quark masses at high energies. At energies 
smaller than the quark masses, the quarks can be inte- 
grated out and the theory looks like the unflavored one 
(a confining SYM-like theory in the present case). This 
suggests, via the holographic radius-energy relation, that 
we can choose, as a model for x(p), just a Heaviside step 
function x(p) = x<d(p — p q ). At the level of the super- 
gravity equations of motion, the only constraint on x(p) 



is that x(p)' > [7]. The exact form of the function 
should follow from the smearing of the massive embed- 
ding associated to the flavor branes. Nevertheless, the 
Heaviside approximation is expected to correctly capture 
the qualitative features of the system under study, as it 
also happens in other similar setups [8]. Since the first 
two equations in (2) are constraints and not differential 
equations, not all the functions can be continuous if x(p) 
is discontinuous. We impose that the metric is continu- 
ous so H, G, a, Y, (p are continuous at p = p q . This is the 
equivalent procedure to matching the couplings when one 
integrates out quarks in the field theory. The only func- 
tion which does not enter the metric is b. In any case, 
what should be continuous - sec cq. (2) - is *fb+ |. The 
solution for b consistent with the differential equations 
and this condition is: 



2p 



sinh(2/9) 



+ X&(p- p q ) 



sinh(2pq ) 
P q ~ P ~ 2cosIi(2p,) 

sinh(2p) 



Below p q , we have the unflavored system, so we can take 
the regular solution, found in Section 8 of the first pa- 
per in [4], depending on a parameter p. Tuning p for a 
given pair (x, p q ), it is possible to get a solution with 
the expected UV which is uniquely determined up to the 
additive constant 4>ir- The system can be solved numer- 
ically. The plots of figure 1 depict a typical case for the 
functions Y, (p. 




FIG. 1: The Y, <p functions in the case x = 0.8, p q = 0.5. 

3. The quark-antiquark potential.- We focus on 
the interaction energy between two test quarks Q, Q put 
as external probes of our theory. The string dual of the 
QQ system is a macroscopic open string attached to a 
probe flavor brane which is very far from the bottom 
of the space (pQ 3> p q ), so that the associated quark 
mass Mq — > oo [9]. The string bends in the bulk and 
reaches a minimal radial position pq. The Minkowski 
separation L between the test quarks, as well as the total 
energy of the system, depends on po- From this we can 
deduce the V(L) relation, where V(L) is the potential 
(renormalized energy), i.e. the total energy minus the 
total quark mass 2Mq. The open string embedding is 
chosen as t — r, y = a, p = p(y) where y G [— L/2,L/2] 
is one of the spatial Minkowski directions. The string 
worldshcct action reads 
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dtdyJg tt (p)[g yy (p) + {d y p) 2 g pp (p)}. 
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Defining / = ^g t t9yy = ^ , g = ^/mg~p~p = 
2 e 0-0/Hy / y v /g s N c a', the (constituent) quark mass is 
given by the energy of a string extended from the bottom 
of the space at p = to the bottom of the corresponding 
flavor branc, m q = J p9 <? dp, M Q = ^7 J PQ g dp. 
We write the length L and potential V as functions of 
Po and with UV cutoff pq as (the subindex means that 
the quantity is evaluated at p = po): 

L(p ) = 2 [ PQ 9fo dp , 

2 r f PQ af f PQ i 

"<»>- w[£ 77%*-/ 

It is not difficult to check the relation 4^ = 4^. 

apo 2na dpi;) 

This for instance means that 4^- and have the same 

dpo dpo 

sign. In some cases they can change sign simultaneously 
at some value of po so the V(L) plot turns around. 

We can now plug our numeric solutions into the def- 
initions of the quark- antiquark potential V(po) and dis- 
tance L(po), to see how they depend on the dynamical 
massive flavors. In doing so we consider a finite (large) 
UV cutoff, and so a finite Mq ^> m q , which we keep 
fixed as m q and x are varied (we will show the results 
for Mq = 200y / g s N c /a' but their qualitative features do 
not depend on this choice). 

4. A first order transition in the potential.- In 
the m q = case, the QQ string breaks at a certain L max , 
where the V(L) plot turns around [4]. At m q ^ 0, this 
turn-around (that one would expect at least in the small 
mass limit), should be followed by a second turn-around, 
such that the potential comes back to an increasing lin- 
ear behavior at large L (i.e. in the far IR, where the 
quarks are integrated out and the theory is SYM-like). 
Our numerical analysis confirms this expectation, adding 
a crucial piece of information: the existence of a critical 
value m c for the dynamical quark mass. For m q > m c , 
the V(L) plot is monotonic, a Coulomb- like behavior at 
small L smoothly joining with the linear one for large 
L. Instead, in the small mass regime, m q < m c , we 
find a first order [13] transition between the Coulomb-like 
and the linear phase, with a corresponding double turn- 
around behavior of V(L) at intermediate distances, see 
figure 2. In the figure, the energetically favored branch 
of the graph is the one with largest L at a given V (the 
solid line). So, there is a discontinuity in the derivative 
of V(L), with a sudden decrease of the "local string ten- 
sion" , such that at a certain energy scale and length the 
increase of the energy of the flux tube with the length 
has a deceleration. The Coulomb-like regime does not 
connect smoothly with the linear regime [14]. 

The existence of a critical mass m c above which there 
is no phase transition, is evident by analyzing the L(p ) 
curve at different values of m q . For example, for x=0.8, 
one can check from plots as those in figure 3 that below 
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FIG. 2: The potential for small dynamical mass m q — 0.01 
and x = 0.8. In the frame, a zoom of the critical region. 
Dashed (dotted) lines are metastable (unstable); cfr. [10] for 
stability issues. Lengths are in units of ^a'g 3 N c ~ 1/Aqcd, 
energies in units of \fa.'g s N c /a' ~ T/Aqcd, T being the 
string tension. 

L 

4.5r , v 

4.25 \ \. 

3.5 ••' ^^^^^ 

I .' 
3.25- . • 
I • 

3- 

2.75- f 

0.20.40.60.8 1 1.21.4 Po 

FIG. 3: Some representative L(po) curves at x = 0.8 and 
m q = 0.01,0.09,0.2 (left to right). Solid, dashed and dotted 
lines correspond to the regions of figure 2. 

m c <~ 0.09, there is a first order phase transition between 
a large po string and a small po one. Actually these plots 
share some qualitative features with the pressure/ volume 
isothermal curves in a van der Waals gas, where the 
phase transition, occurring below a certain critical tem- 
perature, is between a liquid and a gas phase. It might 
prove useful to model the "connected part" of the phe- 
nomcnological potential between two heavy quarks with 
some function, similar to the Gibbs free energy /pressure 
relation in the van der Waals system, exhibiting a discon- 
tinuity in the first derivative for small dynamical quark 
mass. We checked that the quark distance at which the 
transition occurs is a decreasing function of x and an in- 
creasing function of m q . On the other hand, the value 
of the potential at the transition is a decreasing function 
of both x,m q . Finally, m c — > for x — > 0, where one 
recovers the quenched approximation. 

5. String breaking length.- In the presence of dy- 
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FIG. 4: The string breaking length L a f, as a function of: (left) 
the number of flavors x for m q — 0.05,0.1,0.2 (bottom to 
top) and (right) the mass m q for x = 0,0.3,0.8, 1.2 (top to 
bottom). 

namical quarks there is a screening length, defined as 
the length at which the QQ string will break by produc- 
tion of quark pairs. It is determined by the condition 
that the energy of the Wilson loop equals the energy of 
two mesons composed by one test and one dynamical 
quark. In the case at hand, due to the smearing pro- 
cedure, the lighter of such mesons arc nearly massless. 
In fact, in the string language they are localized at the 
intersection of the test brane and the dynamical brane. 
Their mass should scale as Mq/(g s N c ) 1 with g s N c 3> 1 
[11]. This is parametrically smaller than the typical en- 
ergy Eqq(L) = 2Mq + V(L) computed above. Thus, 
the Wilson loop has always enough energy to decay into 
these mesons and the screening length is not a relevant 
quantity. Nevertheless, again due to the smearing, the 
decay in this channel (and in every channel correspond- 
ing to a specific flavor) is 1/N C suppressed. If we ask 
for a decay rate which is not 1/N C suppressed we have 
to consider a larger Q, Q separation, such that the QQ 
state can decay into a sizable fraction of the Nf types 
of heavy-light mesons. Given such a fraction, one could 
study how the separation L s t,, required for the decay on 
all the included flavors to happen, depends on the pa- 
rameters m q , x. Clearly, the choice of the fraction would 
be arbitrary. 

In order to get an idea of the general trend, we choose 
to define as the length at which the Wilson loop has 
enough energy to decay by production of the particular 
flavor which has precisely the same charge (internal or 
R-symmetry quantum numbers) as the test quarks. The 
string dual to the corresponding specific meson, whose 
mass is EqQ = Mq — m q , is stretched only along the ra- 
dial direction of the geometry (1), from the bottom of 
the test brane to that of the "parallel" dynamical one 
[12]. We call L s b the "string breaking length" . The value 
L s b at which the breaking can occur is given by equat- 
ing the energies of the initial and final configurations: 

e qq( l ) - E qQ+ q Q = V ( L ) + 2m q = °- We are neglecting 
the interactions between the two final mesons because 
their contribution is subleading due to the smearing. 
The results are in figure 4. The string breaking length 



is a decreasing function of the number of flavors x: the 
more flavors there are in the theory, the easiest is the 
decay. Finally, the string breaking length is a decreasing 
function of m q : this is a genuine large coupling effect, and 
may change if one chooses a different definition of EqQ. 
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